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FLOW  PAST  A  PLATE  AT  HIGH  SUPERSONIC  VELOCITIES 

V.  I.  Kholyavko 
(Khar 'kov) 


In  accordance  with  the  analogy  of  an  "explosion 
wave"  the  problem  of  high  supersonic  flow  past 
bodies  reduces  to  studying  the  gas  movement  in 
the  transverse  plane. 

The  problem  of  a  flat  explosion  is  equivalent 
to  symmetrical  flow  past  a  flat  blunt  plate.  The 
solution  of  this  problem  has  a  particularly  simple 
form  if  we  disregard  the  effect  of  the  initial 
pressure  (self -similar  solution  of  the  problem  of 
a  "strong  explosion"  [lj). 

The  problem  of  a  flat  explosion  with  subsequent 
movement  of  the  piston  corresponds  to  flow  past 
a  flat  blunt  plate  at  an  angle  of  attack.  In  this 
case,  even  if  we  disregard  the  initial  pressure  the 
problem  is  significantly  complicated,  and  its 
precise  solution  is  possible  by  numerical  methods. 

The  approximate  solution  of  this  problem  can  be 
carried  out  if  it  is  limited  to  small  angles  of 
attack  when  blunting  of  the  plate  exerts  the 
determining  Influence  on  the  nature  of  the  flow. 
Additional  influence  due  to  rotation  of  the  flow 
on  the  plate  (piston  movement)  is  determined  from 
the  equation  for  a  shock  wave. 

We  note  that  an  approximate  solution  to  the 
analogous  problem  of  flow  past  a  blunt  wedge  was 
obtained  by  G.  G.  Chernyy  [2]  by  an  expansion  of 
the  flow  parameters  Into  series  of  powers  of 
(7  -  l)/(7  +  l),  where  7  is  the  ratio  of  specific 
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heats.  The  solution  of  the  problem  concerning 
flow  past  a  blunt  plate  at  an  angle  of  attack 
by  a  method  of  small  disturbances  has  been 
given  [J], 

1.  We  will  Introduce  an  equation  for  determining  the  position 
of  a  shock  wave.  For  this  we  will  use  Integral  ratios  of  the  law 
of  conservation  of  energy,  momentum,  and  mass,  applying  them  to  a 
volume  of  gas  enclosed  between  the  piston  (surface  of  the  body)  and 
the  shock  wave. 

If  we  assume  that  the  width  of  the  gas  layer  is  equal  to  unity 
and  disregard  the  Initial  pressure,  the  Indicated  ratios  have  the 
form: 

5  (¥  +  T=t) dr=*Ejr\ Po  =  7t\pudr<  Pifl-Jp**  (1.1) 

r#  0  **•  r» 

Here  E  is  the  energy  released  during  the  explosion; 
p  is  the  density; 
p  Is  the  pressure; 

u  is  the  velocity  of  the  transverse  gas  flow; 

7  is  the  ratio  of  specific  heat; 

Tq  and  R  are  the  positions  of  the  piston  and  shock  wave  at  time  _t. 
We  will  indicate  parameters  on  the  piston  (surface  of  the  body)  by 
the  subscript  0,  the  undisturbed  gas  by  the  subscript  1,  and  the 
parameters  immediately  behind  the  shock  wave  by  the  subscript  2.  The 
dot  over  rQ  and  R  represents  differentiation  with  respect  to  jfc . 

The  solution  of  system  (l.l)  can  be  carried  out  by  using  some 
distribution  of  parameters  in  the  transverse  layer  between  the  piston 
and  the  shock  wave  [4], 

Let  us  now  turn  to  the  third  equation  of  (l.l).  It  follows 
from  this  equation  that  regardless  of  what  law  of  density  distribu¬ 
tion  in  the  transverse  layer  p(r)  is  taken,  the  magnitude  of  the 
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integral  on  the  right  side  will  remain  constant.  In  particular, 
this  will  also  be  the  case  if  the  entire  gas  mass  is  concentrated 
about  the  shockwave  [2,4]. 

In  this  limiting  case  the  density  distribution  can  be  represented 
in  the  following  forms 

P  (r)  =  PiA&  (r  —  R) 


where  6(r  -  R)  is  the  Dirac  function  (delta-function). 

After  this  it  is  not  difficult  to  determine  the  following 
integrals  from  System  (1.1) 


r.  r> 


/ 


where  Ug  is  the  velocity  behind  the  shock  wave  [l]. 

Disregarding  the  pressure  change  across  the  layer  [2]  we  will 
also  calculate  the  integral 

If  we  now  substitute  these  approximate  relations  into  the  first 
and  second  equations  of  (1.1),  we  then  obtain  the  pressure  on  the 
piston  (surface  of  the  body) 


Po 


2Pi 

T+l 


(A*  +  RR) 


(1.2) 


the  equation  for  determining  the  position  of  the  shock  wave 

t 

Rfr  +  Ii'-R  -  r 0  {ir-  +  RR)  =  |  +  (T  -  1 )  J  r0(Rt  +  RR)  dt  ( 1 .  J  ) 

0 

In  a  case  of  a  flow  past  a  plane  at  an  angle  of  attack,.  rQ  » 
aVt  (V  is  the  incident-flow  velocityj  a  is  the  angle  of  attack; 
for  the  upper  surface  the  'value  a  is  taken  with  the  minus  sign)  and 
Eq.  (1.5)  assumes  the  following  form: 


-> 
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AVf*  +  -  *v  (it-  +  )t  =  £.  +  (T  _  i)  am/  (1.4) 

2.  We  will  examine  the  particular  cases  of  solving  Eq.  (1.4). 
After  setting  a  *  0  we  obtain  a  solution  of  the  problem  of  a  "strong 
explosion"  (R  =■  RQ)  which  corresponds  to  symmetrical  flow  past  a 
blunt  plate 

+  Ii*Ra  =  'L =1  JL 

Prom  this. 


and. the  pressure  magnitude  according  to  Eq.  (1.2) 


.  r*  (t  + 1)  (r  —  D*  T/»  t l:  S'1'  ,->i. 
Pi  Px  LT  (3T  -1?  J  Vi?rJ 


(2.1) 


(2.2) 


In  order  to  now  obtain  flow  past  a  flat  blunt  plate  it  is  sufficient 
to  substitute  into  Eqs.  (2.1)  and  (2.2) 

E  =  —  Cx^-d,  r  =  -i- 

where  C  Is  the  blunt  drag  coefficient,  d  is  the  thickness  of  the 

A  " 

plate;  as  a  result  we  obtain 


jk.  r 9 

d  ™  [  16 


(T+1)»  (T-l) 


M-f)' 


•(r-  i)(r-))»  T/.  1  \-v. 

-i)5  J  '*'[-) 


(2.5) 

(2.4) 


Formulas  (2.5)  and  (2.4)  coincide  with  the  solutions  obtained 
by  G.  G.  Chernyy  [2j.  This  same  work  gives  a  comparison  with  the 
calculations  by  a  precise  theory  of  a  "strong  explosion." 

We  can  show  that  (2.1)  is  the  solution  of  Eq.  (1.4)  not  only 
when  a  -  0  but  as  well  as  for  any  other  a  /  0  if  the  limiting  case 
t -4  0  is  considered. 
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In  another  limiting  case  t-*  »  the  solution  of  Eq.  (1.4)  and  the 
pressure  magnitude  have  the  following  farms,  respectively: 


D  T  *M  i  T — '  1  E 
R  —*vt+  Pl 


^r=(r  +  1)«* 


(2.5) 


The  pressure  magnitude  determined  by  (2.5)  coincides  with  the 


exact  solution  for  a  tapered  plate. 

Thus,  if  close  to  the  leading  edge  of  the  blunt  plate  (t  -  x/V-» 
-»  0)  the  flow  is  determined  mainly  by  blunting,  then  far  downstream 
(t-»  eo)  the  effect  of  blunting  diminishes  and  the  pressure  magnitude 


tends  toward  the  value  occurring  at  a  tapered  plate. 

In  order  to  obtain  a  solution  to  (1.4)  for  small  but  finite  t, 
we  assume  that  R  -  R0  +  R^,  where  R2  is  defined  (2.1)  and  R^  is  the 
addition  due  to  rotation  of  the  flow  (piston  movement).  Here  no 


assumptions  are  made  relative  to  the  magnitude  of  R^;  however  the 
following  conditions  must  be  satisfied:  l)  when  t  ->  0,  R^-*  0  and 
the  pressure  magnitude  should  tend  to  a  value  (2.2);  2)  the  appear¬ 


ance  of  additional  magnitude  R^  is  connected  only  with  the  presence 
of  an  angle  of  attack  (piston  movement);  in  other  words,  parameter 
a  should  definitely  enter  into  the  solution  of  R^. 

If  we  now  assume  R^ —  tm  and  set  R  -  RQ  +  R^  into  Eqs.  (1.2) 
and  (1.4),  we  find  that  the  first  condition  will  be  fulfilled  when 
m  >  2/5,  and  the  second  when  m  >  5/6. 

Leaving  in  converted  equation  (1.4),  i.e. ,  after  substitution 


of  the  relation  R  =  Rq  +  R^.  into  it,  the  terms  containing  t_  to 
the  lowest  power  (when  m  >  5/6  these  will  be  the  terms  with  t1^ 
and  tm-2/,,^),we  obtain  an  approximate  equation  for  determining  R^ 


l'Jtl  ~  T  r+T  Uii  +  T  f+T  Rl  m  T  (3T~ 


(2.6) 
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The  equation  under  consideration  is  an  Euler  equation  with  the 
right  side  [5]. 

The  general  solution  of  this  equation  has  the  following  form: 


+cjf+  +  Itu 

where  c Cg  are  arbitrary  constants,  m^,  mg  are  roots  of  the  charac¬ 
teristic  equation,,  and  is  the  particular  solution  of  (2.6). 

Solution  of  Eq.  (2.6)  leads  to  negative  values  of  nu  and  r:.0 
(7  >  1)  and,  therefore,  in  order  to  satisfy  the  boundary  conditions 
it  is  necessary  to  assume  that  c^  =  Cg  -  0. 

Thus,  only  the  particular  solution  of  (2.6)  remains.  It  is 
not  difficult  to  see  that  it  will  be  as  follows: 


n  _(r+D(3r-2) 

■  1 - iifr-i) 

Finally,  the  solution  of  Eq.  (1.4)  for  small  values  of  jt  has 
the  form 


*-Kjug$=flr(*r*+ 


(T  -f  1)  (37  —  2) 
2(7T-i) 


<xVt 


(2.7) 


By  substituting 


±c  liHd 
2  c*  2  a  * 


Formula  (2.7)  Is  reduced  to  form 


is(il)  =■  li(x))v.  +  *(.) 


s  (*)  =  _ 


(2.8) 


d 


It  is  completely  obvious  that  the  obtained  solution  of  (2.7)  or 
(2.8)  is  correct  not  only  during  flow  past  a  blunt  plate  at  an  angle 
of  attack,  but  also,  e.g.,  during  symmetrical  flow  past  a  blunt  wedge. 

Figure  1  shows  the  comparison  of  results  of  calculation  by 
Formula  (2.8)  (dashed  line)  with  the  solution  of  Q.  0.  Chernyy  (7  -  1.4) 
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for  a  case  of  flow  past  a  blunt  wedge  [2]  (on  the  ordinate  axis 
X(R)  -  2z(R). 

Using  (2.7),  we  determine  by  Formula  (1.2)  the  pressure  on  the 
surface  of  the  blunt  body 

Tp* 

-A  (?)*-*'•  +  ft  (y)*-1'*  +  ft  (Y)  (2.9) 


Here 


/» (r) 


•(T-HHr—  D'T/. 

.  36(3t-1)>  ’ 


/*  (r) 


A(r) 


(r  +  i)(3r-2)« 
2(7y-1)» 


5  3r — 2  rcr  -f- 1  )*  f  r — i )  t/» 

3  7t  —  1  L  8  (3f  —  1) 


Fig.  1. 


Fig.  2. 


In  particular,  for  7-1.4  the  following  result  is  obtained; 

P’  “  p^pr  -  O.lz-V.  +  0.206*-%  +  0.075  (  2. 10) 


For  comparison  let  us  Introduce  the  formula  which  was  obtained 
by  the  method  of  small  disturbances; 

p'  -  -^3-  -  0.085*-%  +  0.194*-%  (2.11) 


Figure  2  shows  the  results  of  calculations  by  Formulas  (2.10) 
(dashed  line)  and  (2.11)  (solid  line).  The  small  circles  and  the 
small  dashes  indicate  Bertram's  data  which  he  obtained  by  the  method 
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c.  characteristics  for  a  -  5°,  10°,  with  -  20  [?]. 

3.  Let  uq  tprn  to  determination  of  the  aerodynamic  characteris¬ 
tics  of  a  flat  blunted  plate  with  flow-past  at  high  supersonic  veloci¬ 
ties  at  small  angles  of  attack  a. 

In  accord  with  Formula  (2.9)  the  coefficients  of  relative  pc 
at  the  lower  and  upper  surfaces  are  determined  by  the  following  ratios 


-'A (?) ^ (tP  +  V, (T) c,v. fif'a  +  2/, (r) a2  (3.1) 

P- -  -  2A  (T) -  2A (r) *«•'•  (fp*  +  2/a  (r) oc»  (?<2) 

From  (3.2)  it  follows  that  on  the  upper  surface  there  y 

exists  such  a  point  0,  at  which  the  pressure  will  be  equal  to  zero. 

From  physical  considerations  this  point  can  be  examined  as  a 
point  in  which  breakaway  occurs.  If  we  disregard  the  third  term  as 
a  small  magnitude  compared  with  the  other  terms,  from  (3.2)  we  obtain 


9  (T-OffT-l)*  c* 
250  (T  + 1)  (»7  -  1)  (»7  -  2)*  a* 


(3.3) 


Fig.  3. 


Fig.  4. 

Figure  3  shows  the  dependence 
(3.3)  when  7  -  1.4  and  Cx  ■  1.84 
(flat  nose).  As  is  evident  from 
Fig.  3,  with  an  increase  of  the 
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angle  of  attack  the  breakaway  point  moves  forward  toward  the  leading 
edge  of  the  plate.  At  large  angles  of  attack  the  breakaway  point  is 
located  on  the  leadings  edge,  and  the  pressure  on  the  whole  upper 
surface  will  be  equal  to  zero.  A  similar  phenomenon  occurs  during 
flow  past  a  plate  according  to  the  Newton  system. 

If  we  limit  ourselves  to  an  examination  of  such  angles  of  attack 
for  which  point  0  is  positioned  behind  the  plate,  then  the  lift  coef¬ 
ficient  is  defined  as: 

1 

*  =  | ~  p-]  d  (7)  =  6A  (r)  cx'u 

Here  la  is  chord  and  <1  is  the  thickness  of  the  plate. 

We  obtain  the  total  lift  coefficient  if  the  pressure  on  the 
leading  part  of  the  plate  is  taken  into  consideration: 

(3.4) 

Figure  4  shows  the  dependence  (3.4)  constructed  for  7-1.4  and 
Cx  -  1.84  (curve  l).  For  comparison,  we  have:  the  lift  coefficient 
of  a  blunt  plate  obtained  without  taking  into  account  the  effect  of 
blunting  on  flow  past  a  plate  (curve  2): 

CI/=(Tt1)»’  —  Cx«  y 

the  left  coefficient  of  a  tapered  plate  (curve  3) 

c,,  =  (r-M)*3 

As  was  to  be  expected,  at  large  angles  of  attack  the  effect 
due  to  bluntirg  diminishes  and  the  characteristic  c_  -  f(a)  of  a  blunt 
plate  approximates  those  values  obtained  for  a  tapered  plate. 

The  coordinate  of  the-  center  of  pressure  X  on  the  flat  blunt  plate 
relative  to  the  leading  edge  is  defined  by  the  following  expression 
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X 

V 


(5.5) 


(J  (/>*  -  P-)  j  d  t)  ($  (/>♦  -  P- )  *t)  “  0,4 

o  u 

Thus,  due  to  blunting  the  pressure  center  of  a  plate  shifts 
forward.  With  an  increase  of  the  angles  of  attack  the  pressure  center 
will  move  into  a  position  corresponding  to  a  tapered  plate. 
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